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Advances in material growth methods have renewed the interest in localization of one-dimensional 
systems in the presence of scale-free long-range correlated disorder potentials. We analyze the 
validity of single parameter scaling for the /3 function away from the band center, in the presence 
of correlations. A renormalized disorder strength emerges reducing the regime of validity of the 
single parameter scaling hypothesis. Analysis of localization lengths and participation ratios leads 
to correlation dependent critical and fractal exponents, consistent with the extended Harris criterion. 



The Anderson model of disordered one-dimensional 
(Id) systems provides valuable insight into the elec- 
tronic transport properties of a wide variety of mate- 
rials such as DNA[TH3], polymers 0HS], random media[7] 
and nanowires[5], for which the scaling theory of local- 
ization assumes exponentially localized eigenstates with 
vanishing conductance [5]. Historically, the successful ap- 
plication of scaling theory to systems with short-range 
correlations[10 was followed by the development of cor- 
related disorder mo dels [IT]. These models fall into three 
basic categories: those with quasi-periodic real space 
order [T2l 113) . random disorder amplitudes chosen from 
a discrete set of values [TH [15], and models with spe- 
cific long-range correlations jTS], In the first two cate- 
gories, a discrete number of eigenstates are predicted to 
remain extended. The last class of models however, pro- 
duce a mobility edge that separates regions of extended 
and localized states in the Hilbert space, rendering a 
metal-insulator transition in one-dimension. As a conse- 
quence, much work has been devoted in the last few years 
to understand the effects of different long-range correla- 
tions and, in particular, those with power-law spectral 
functions, generated by self-affine potentials[T7]. Disor- 
der potentials of this class however are not easily re- 
alizable in controlled experimental setups, thus several 
authors have focused on potentials with power-law cor- 
relations in real-space. These have already appeared in 
experiments in ultra-cold atom systems [T51 119) . where lo- 
calization physics was achieved using speckle potentials, 
known to possess long-range correlations. Furthermore, 
advances in materials growth methods such as molecular 
self-assembly on surfaces or wires [23 [2T| , and the produc- 
tion of rippled graphene flakes and ribbons [221124) . makes 
it possible to obtain systems where real-space correlated 
disorder potentials can be manipulated, making the un- 
derstanding of these models of experimental relevance. 

Transfer matrix studies for real-space power-law cor- 
related models have shown, for all disorder strengths: i) 
complete localization for all values of the power-law ex- 
ponent a, ii) anomalous enhancement of localization for 



band edge states, iii) more extended band center states as 
compared to those in the presence of short-range corre- 
lated potentials 25, [26]. Furthermore, a scaling analysis 
of a weak disorder expansion predicts a rescaling of the 
disorder strength W — > W e ff that takes place near the 
band center, where W eff = ^/ S(2k)W\T7\, and S(2k) 
is the 2k component of the spectral density associated 
to the correlated potential. Finally, it was shown that a 
white-noise disorder model reproduces the effect on the 
band edges 26J, suggesting a different scaling behavior 
of the potential amplitudes in this energy region. Note 
that these studies have not addressed the anomalies at 
the band edge and centers found in studies with uncor- 
rected potentials [28H30] 

Although it is well established by now that all eigen- 
states are localized for these correlations, other effects 
with direct relevance for applications have not been fully 
addressed up to date. In this work, we focus on three 
of them: i) validity of single parameter scaling (SPS), ii) 
crossover between the two predicted scaling regimes, and 
iii) nature of localized states as described by the local- 
ization length critical exponent v and the participation 
ratio fractal exponent D. 

First, we analyze the single parameter scaling hypoth- 
esis for the j3 function (an issue that remains contro- 
versial at present [T7j ) . Second, we calculate the local- 
ization length £ and study the crossover that separates 
eigenstates into two regions: those affected by W e ff and 
those affected by an effective white-noise disorder. We 
characterize it by analyzing the change in v. Finally we 
calculate the inverse participation ratio and obtain the 
corresponding values for D for localized wave functions. 
Our results suggest that v and D acquire a dependence on 
the correlated disorder power-law exponent a, consistent 
with predictions by the extended Harris criterion |31U32| . 

Model: We begin with the Anderson hamiltonian for 
a one-dimensional chain of N = L/a sites in real space, 
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with L the length of the chain and a the lattice constant: 
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Here e n is the random on-site energy , t is the hopping 
energy, and (c n ) is the fermion creation (annihilation) 
operator at site n. e n is obtained from a gaussian distri- 
bution satisfying the conditions (e„) = 0, (e 2 /) = W, and 
with a normalized correlation function 
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where n is the distance between sites, and a determines 
the correlation strength. Although several choices render 
the same asymptotic form l/n" as n — > oo|26|. short 
range details are irrelevant in the thermodynamic limit 
which are the focus of this work. 

Methods: To generate a stochastic set of real on-site 
energies, a set of random complex numbers in reciprocal 
space (k-space) is generated, with the condition ik = £*_ k 
(e~fc is the complex Fourier component of the real space en- 
ergy e n ). These variables are chosen from a gaussian dis- 
tribution with a standard deviation equal to the spectral 
density S(k) = (|efc| 2 ). A Fourier transform of Eq.([2]), 
renders S(k) in a computationally friendly form: 
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Here, FQ and Re() denote the Fourier transform and real 
part of the argument respectively. The discrete spectral 
density is obtained with a Fast Fourier Transform (FFT) 
algorithm, a random set {e~k} is generated, and the set of 
{e n } is obtained via the inverse Fourier transform. 

Scaling function: The validity of the SPS hypothesis 
for the dimensionless conductance of systems with long- 
range correlated potentials is still an open quest ion |17j. 
For example, for the disorder potential considered in 
Ref.[I5], the exponent determines an effective length 
scale beyond which the real-space correlation function be- 
comes negative[33,,[34]. To answer the question, for l/n a 
correlated potentials, we use the relation Var(ln g) — 
— 2(lng) between the first (average) and second (vari- 
ance) cumulants of the "typical" conductance lng, found 
to hold for uncorrclated disorder potentials 35, 35J, as the 
criterion for the validity of SPS. 

The dimensionless conductance g is equal to the trans- 
mission function T obtained from advanced and retarded 
Green's functions by T = Tr(T L G r c T R G a c ) where: 



G r c ' a (E) = 
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Here E is the energy, Hd the hamiltonian of the disor- 
dered system, and E2°r se ^ ener gi es due to left and right 
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FIG. 1. Scaling of the average value of the typical conduc- 
tance and it's relation to the variance. Data is averaged over 
1000 disorder configurations. Symbol size represents statisti- 
cal error, (a) Raw conductance data for L/a — 2 9 -2 17 and 
various parameters, (b) Data collapse by rescaling of L by 
£ as described in the text. Notice the universal slope for all 
values in a.) (c) Relationship between variance and average 
for L/a=2 17 . Dashed black line shows Var(lngr) = 2(— lngr). 



leads. The hybridization functions Tl and Tr are defined 
byr^^E^-E^]. 

Usually Ej'fl is found numerically through decimation 
procedures |37j. however for semi-infinite one-dimensional 
leads, an analytic expression can be derived [55]: 
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To calculate the inverse of G r c ' a {E) we use recursive meth- 
ods that require an order ~ N of operations [39] , allowing 
to maximize the number of disorder realizations and sys- 
tem sizes. Results from these calculations are plotted in 
Fig. [l] Panel (a) and (b) show the collapse of each curve 
after a rescaling of L by £ (with £ = limL^oo ) , in- 
dicating d }\^ 9 ]] = (In 9) ■ The value of the /3 function is 
negative confirming localization physics. Panel (c) shows 
the relation between the first two cumulants. The data 
suggests SPS is valid for low energies but violated near 



the band edge for ^ < 1 and for all energies if 4— > 1. 

We define E sps as the energy where SPS is first vio- 
lated, for each value of a and W/t. Fig. [2] shows the 
dependence of E sps on a (a) and W/t (b). By defining a 
rescaled variable all curves collapse as shown in (c). 
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FIG. 2. Energy E sp3 separating valid (E < E sps ) and vio- 
lated (E > E sps ) SPS regions as function of a (panel (a)) 
and disorder strength (panel (b)). Panel (c) shows E aps as 
function of the rescaled variable The black dashed line 
represents the fit of Eq. rol 



Numerical fits of the data lead to the expression 
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reminiscent of the relation found in Ref. [3D] for the energy 
at which the perturbation expansion of £,{E) done with 
uncorrelated potentials near E ~ 2 ceases to be valid pfO"]. 

Localization Length: Next, we apply transfer matrix 
techniques to obtain £ as a function of energy [25l [26] 
by calculating the lowest Lyapunov exponent from the 
successive multiplication of a position dependent transfer 
matrix 
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We show a typical data set in Fig. [3]ja). The values plot- 
ted are obtained by calculating £ for different systems 
sizes and extrapolating to L — > oo. As predicted from 
scaling arguments, and partially confirmed in transfer 
matrix calculations [25J [27J , two regions are distinguished 
with enhanced and suppressed localization. These re- 
gions are separated at the crossing energy E c /t ~ y/2, 
defined as the energy where £ is independent of a for a 
fixed value of W/t. These results indicate that correla- 
tions produce more extended states in general (i.e. the 
region of 'band-center' states extends closer to the band- 
edge) with localized states 'squeezed' into the band edges. 
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FIG. 3. (a) Localization length vs energy for system sizes 
L/a — 2 11 -2 19 and 1000 disorder configurations. Eigenstates 
become more extended with decreasing values of a in region 
I. At E c /t ~ V2 (dashed line), £ is independent of a. In 
region II, there is enhanced localization. (Positive values of 
E are plotted only), (b) Dependence of crossing energy E c /t 
vs disorder strength W/t. Solid, colored region corresponds 
to energies with enhanced localization while dashed region 
corresponds to more extended-like states, (c) Localization 
length critical exponent v as a function of a for W/t = 1 
showing different functional dependence for a < and > 1. 
Symbol sizes represent statistical error bars. 



They also suggest a continuos enhancement of extended 
states as a is reduced with the singular limit of perfect 
extended states for a = 0. The energy E c , that depends 
on the disorder strength as shown in Fig. [3^b), separates 
the region of validity for the two scaling regimes proposed 
in Refs. [2S1[27]. A quantitative analysis of the critical 
exponent v obtained from £ ~ \E — E c \~ v (with E c = 
in this model), produces the results shown in Fig. |3jc) 
for W/t = 1. The data suggests v = v(a) with a func- 
tional dependence that changes at a = 1, consistent with 
the extended Harris criterion for classical [31] and quan- 
tum systems [32] . Similar behavior is obtained for other 
values of W/ 1 not shown. 

Participation ratio: To get insight into the nature of 
localized eigenstates we performed an analysis of the par- 
ticipation ratio. This quantityjH] provides information 
about real space fluctuations of eigenstate amplitudes. 
For a state m it is defined by: 



Am — 



J2n=o KMn)| : 



N E^hMn)! 4 



(8) 



where tp m (n) denotes the amplitude of the eigenstate at 
position n. If the state is localized, A m approaches zero 
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FIG. 4. Participation ratio A m for states at energies E/t = 0.1 
(a) and 1.7 (b) for L/a = 100-1000 averaged over 100 disorder 
realizations (error bars obtained from disorder averaging) and 
W/t=l. Largest value for £/a ~ 30. Solid lines are fitting 
functions of the form A = C(a/L) D . Error bars are obtained 
by error propagation from proposed fitting function. Panels 
(c) and (d) show values of fractal exponent D for different 
values of a. A change at a = 1 is observed, consistent with 
the SPS results. 



as A m rgr^vD as the system size increases. Here D is 
the fractal exponent, that characterizes the fragmented 
nature of localized eigenstates [4"2HH] . 

An exact diagonalization of Eq. ([I]) renders the eigen- 
states and the numerical evaluation of A m as a function of 
system size. Fig. [4] shows results for two different states, 
chosen near the band center (a) and the band edge (b) 
respectively for W/t — 1 as representative data. Con- 
tinuous lines correspond to a fitting function c (L/a)~ D 
with D and c a fitting parameters dependent on W/t, E/t 
and a. In panels (c) and (d), the fractal dimension D 
is plotted versus the correlation length exponent a. The 
data suggests the existence of two regimes for a > 1 and 
a < 1. In the first regime D is independent of a, a 
result expected for uncorrelated disorder potentials |4"5] . 
However, for a < 1, D acquires an a-dependence that ap- 
pears stronger for band-center states than for band-edge 
states. The decrease of D with a indicates more extended 
eigenstates (less fragmented regions with non-zero wave 
function amplitude) near the band center. The partic- 
ular dependence of D with a is also reminiscent of the 
relation found between v with a in the previous section. 

Conclusions: We have shown that in the presence of 
real space scale-free correlated disorder potentials with 
power-law (positive) correlations, the SPS hypothesis is 
violated for energies between the band edge and center, 
when the ratio Y~ < 1. For j- > 1, single parameter 
scaling is violated for all energies. The a-dependent en- 
ergy threshold at which SPS is violated, E sps , is consis- 
tent with the value for which the effects of the band-edge 



anomaly appear in the expansion for £ in terms of W/t, 
thus linking the SPS violation to the non-analyticity of 
£. While this could suggest that correlated potentials 
are equivalent to uncorrelated ones with a rescaled disor- 
der strength W' ~ W/a, the anomalous increase of the 
localization length for states with energies E < E c is in- 
consistent with this picture (with E c the crossover energy 
between two scaling regimes). In fact, correlations affect 
the nature of localized states depending on their energy, 
as shown by the variation of v and D as functions of 
a. In particular, they produce an increased number of 
extended-like states for E < E c , more pronounced when 
a < 1, indicating a continuous interpolation to perfectly 
extended states in the limit a = 0. The clear departure 
from SPS with the appropriate degree of correlations, 
ubiquitous in a variety of real systems, calls for the need 
of careful interpretation of experimental data in systems 
with these types of disorder realizations. 
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